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Abstract 

In this paper we investigated the Drell-Yan process of the light neu- 
tralino pair = 1,2) productions at hadron colliders. We 
studied the dependence of the coupling properties of two light neu- 
tralino on the three SUSY Lagrangian parameters Mi, M2 and /1, 
and find that the production rate of X1X2 P arr be dominated by 
the Higgsino-like couplings under the condition /1 <C M 2 and \i ~ Mi, 
while the productions of X1X1 an d X2X2 pairs will be enhanced by the 
gaugino-like couplings under Mj <C /1. For the Higgsino-like Xi 2 n eu- 
tralinos, the cross section of X1X2 production at the LHC can reach 
100 fb. 
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1. Introduction 

With various theoretically appealings, the Minimal Supersymmetric Standard Model 
(MSSM) is considered as the most potential choice for new physics beyond the Standard 
Model(SM). However, due to the absence of supersymmetric particle discovery in the en- 
ergy range of current colliders, the supersymmetry (SUSY) is clearly not an exact symmetry 
of nature, and therefore must be broken. If R-parity conservation is hold and the SUSY 
breaking is around TeV scale, the first two neutralinos Xi,2 are expected to be light, and the 
lightest one Xi is probably the lightest supersymmetric particle(LSP) in the SUSY models. 
Neutralinos would be the most promising particles for the first experimental SUSY test. Due 
to the high luminosity and large energy, the Large Hadron Collider(LHC) has the ability to 
produce neutralino pairs xlxfj- 111 proton-proton collisions, XiX'j pairs can be produced via 
Drell-Yan subprocess and gluon-gluon fusion. Although the anti-quark luminosity in distri- 
bution function of proton is much lower than gluon, yet the cross sections of the neutralino 
pair productions via Drell-Yan mechanism are competitive with those from the gluon-gluon 
fusion, since the former mechanism of neutralino pair productions are accessible at the tree 
level. These facts make the production rates in Drell-Yan process are competitive with or 
even larger than those in gluon-gluon fusions. Therefore, it is significant in theoretical and 
experimental studies of XiX^j P a h productions via qq collisions at LHC. 

There are four SUSY model parameters involved in chargino/neutralino(ino) sector, i.e. 
gaugino mass parameter Mi and M 2 , Higgs mass parameter /i and the ratio of vacuum 
expectation values(VEV) tan/3 = v 2 /vi. These four parameters are essential ingredients 

1 



of the model: the two gaugino mass terms quantify the supersymmetry soft breaking of 
SU(2) x 17(1) subgroup, and together with fi and tan/? determine all the phenomenology 
properties in the ino sector. In order to extract the ino physical quantities (i.e. masses, 
mixings and couplings) from experimental observables (such as cross sections, branching 
ratios and L-R asymmetries), one needs to study thoroughly the relations of these quan- 
tities and observables with the fundamental SUSY parameters. The discussion of the pair 
productions of neutral gaugino and Higgsino current eigenstates can be found in Ref. |T|. 
However, the discussions were not based on physical neutralino particles, and thus could 
not really provide the sufficient information for the investigation of SUSY particles. In this 
work, we studied three typical pair production processes of neutralino mass eigenstates in 
proton-proton collisions, and investigated the dependence of their cross sections on SUSY 
Lagrangian parameters. The paper is organized as follows: in the section below, we present 
some general discussions on the diagonalization of ino mass matrixes. Sec. Ill gives the 
analytical tree-level formula of the cross sections of subprocesses qq — > XiX*j- m Sec. IV, we 
study the production rates of Xi X°j pairs in both the subprocesses and the parent processes, 
and discuss the dependence of the production rates on the SUSY model parameters. Finally, 
a brief summary is presented. 

2. The MSSM parameters in chargino/neutralino sector 

The charginos xt (j ; = 1> 2) mass matrix in the current eigenstate basis has the form 




(2.1a) 
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and the neutralino x° (2 = 1 , 4) mass matrix is 

/ Mie^ 1 — mz sin Ow cos (3 mz sin 6w sin (3 \ 

y _ M 2 m z cos 6 W cos (5 —m z cos sin P 

— m z sin 9 W cos /3 cos^ cos/5 |/i|e l< ^ 

\ sin 9 W sin /3 — m z cos 9 W sin f3 |/i|e l< ^ / 

(2.16) 

Without loss of any generality, \i and Mi can be complex which can introduce CP-violation, 
and M 2 would be set to real and positive for its phase angle can be rotated away by the 
field redefinition and thus absorbed into </>^ and 4>i- The diagonalizations of matrix X and 
Y, which give all the ino masses and their couplings, can be achieved with 

U*XV^ = X D 

N*YN^ = Yd (2.2) 
In this work, we consider the ino sector with the following assumptions: 

• For simplification, CP-conservation is hold, namely 0^ = 0i = 0. 

• The physical signs among Mi, M 2 and /j, are relative, which can be absorbed into 
phases M and <f>i by redefinition of fields. Thus, Mi, M 2 and \i are chosen to be real 
and positive, i.e., Mi,M 2 ,/i > 0. 

• tan/3 is always set as an input parameter, supposed that adequate information about 
tan p could be obtained from some other experiments beyond ino sector. 

With the above assumptions, generally we have two ways to choose the input parameters in 
our calculation. One can employ the scenario of taking M 1 , M 2 , \x and tan/3 as input pa- 
rameters, and then get all the physical ino masses and the matrix elements of U, V and iV as 
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outputs. The detail analyses can be found in Ref. || Since the Lagrangian parameters 

should be extracted from the physical quantities, one can also choose an alternative way to 
diagonalize the mass matrix Y by taking any three physical ino masses together with tan (5 
as inputs. This strategy will provide the other three ino masses and all the mixings and 
couplings as outputs. There are several scenarios about the choices of three ino masses 0. 
On the consideration that adequate information from chargino sector will give out the two 
mass values of xf 2 , and from the energy distribution of the final particles in the decay of X12 
at least one mass value of m^o can be measured, in this paper we take two chargino masses 
m-+ and one of the neutralino masses m^o as inputs. In this way, the two fundamental 
SUSY parameter M 2 and fi can be figured out from the chargino masses by using following 
formula 

(M 2 ,fj,) = M± = X - (^(m-+ + m~+)2 - 2m 2 w (l - sin 2/3) ± yj(m^ - m^f - 2m 2 w (l + sin 2/3)) . 

(2.3) 

The ambiguity in determining fi and M 2 from the two-fold values of M± can be cleaned 
up with some favourable measurements on the chargino phenomenology. Practically, if the 
behaviours of charginos and their couplings are Higgsino-like, there exists M 2 ^> fi, while 
if they are gaugino-like, we have M 2 <C fi. As to Mi, it is usually a free parameter in 
the MSSM, if we have no further assumption. In the second input strategy, Mi can be 
determined by one of neutralino masses m^o alone, when M 2 and /1 are given by Eq.(2.3). 

Mi = - (t2?tm (m~o - /i 2 )(m ? o =p M 2 ) ± 2m 2 ~ Q m 2 z - 2m^M 2 m 2 z sin 2 9 W + 



4 



/i(2m^o =f M 2 )m| sin 2(3 ± //M 2 m| sin 2/3 cos 26^} / 
|(/x 2 - m|o)(m^o =F M 2 ) + m| cos 2 ^(m^o ± /j. sin 2/3)| 

(2.4) 

Generally there are two solutions for M 1; and one is positive while the other is negative. 
With the assumptions mentioned above, we take the positive value solution from Eq.(2.4), 
and fix Mi definitely. Therefore, from the input mass values of two charginos ("\+ 2 ) and 
one of the neutralino mass(m x o), one can extract all the SUSY mass parameters Mi and 
M± which denote M 2 and \i alternatively, and then figure out the mass spectra and all the 
couplings of neutralino sector consequently. In the following calculation and discussion, we 
will adopt both input strategies discussed above. 

3. Neutralino pair productions in qq collisions 

The Neutralino pair productions processes via the collisions of quark and anti-quark in 
protons, can be expressed as 

q(pi)q(p2) - x-ih^ih) = 1,2,3,4) (3.2) 

where p\ and p 2 represent the momenta of the incoming quark and anti-quark, and k\ and 
k 2 denote the momenta of the two final state neutralinos, respectively. The Mandelstam 
variables s, i and u are defined as s = (k\ + k 2 ) 2 , i = (pi — &i) 2 , u = (pi — k 2 ) 2 - The relevant 
Feynmann diagrams are drawn in Fig.l. The interaction Lagrangians involved are listed as 
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§ E 

COS Uyy 
Z COS o\y 

Cqqx° = -Vig ■ q[aqki p L + a^ ki P^xl% ( 3 - la ) 

with the coupling constants 

a| g = Q g sin 2 ^, <=~(-lf« 3+1 / 2 

= iV i|4 iV* 4 - iV ii3 iV*3 
%ki = 2mw(3 g R * q k ' 1 ~ ^QwQqNliR* fej2 
af w = (T g 3 iV li2 - tan^(T g 3 - Q q )N hl )R* q M + k,2 (3.16) 

Here Pr^ = (1 ± 7s)/2. i? g denote the squark transformation matrix, q = 1,2 denote 
up-type and down-type quarks respectively, and 

P q = \ Sin ^ g = J (3.1c) 
^ y 1 cos p q = 2 v y 

The corresponding Lorentz invariant matrix element for the tree-level process is written as 

M = M s + M { + Mu (3.3a) 

where 

M- s = %9 \ u(h)^[O l iP L - OfP R ]v(k 2 ) ■ v(p 2 hA4q + 4qP L H Pl ) 

cos "vv " ) 
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t ~ ml 



«(fci)(a^P fl + a%P L )u( Pl ) ■ v(p 2 )(a L qkj P L + a« kj P R )v(k 2 ) (k = 1, 2) 



A4« = (-l)* 



u(fc 2 )(aJ£P fl + aJgP L )«(pi) • v(p 2 ){a L qli P L + a^P il )^(A; 1 ) (/ = 1, 2) 

(3.36) 



Here we take simply vanished masses of the first generation quarks, i.e. m u = ma = 0. The 
relative sign (—1)^ of M u towards Aif and M$ is merely due to Pauli statistics. 
The corresponding differential cross section at tree level can be written as 



da _ IN 1 Sii g% 
dQ 49T 32vr 2 s 2 



(hs + I tt + hu + 2L si + 2{-l) 5 -h u + 2(-l) s -I iu ) 



(3.4a) 



where the factors |, ^ and (\) 5ij are the initial spin-average, color-average and final 
identical-particle factors respectively. The squares of matrix element have the form as 



,S 2 



cos 4 9w(s — m z 



T~^\"Z< I + "/.;"/.,, + 2 a Zq ) 



O'iOT ■ (« - i)(m% - i) + (ml - u)(m% - u)) - (0% 2 + Of 2 ) • m^s 



hi = 



(t - m\)(t - m 2 qy ) 



(a L q i>Atj + a fi'j a S)( a iu a ivi + «*J • (m\ - u)(m\. - u) 



l st 



_9 



cos 2 6 w (s - m 2 z )(t - m 2 -J 



'AA^P'i - (4 q + a L z q Kki*%jOr) ■ (ml - t){m% - t 



(( a Zq + a Zq) a fki a fkjOz a Zq a qki a qkjOz ) ' m Xi m Xj^ 
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((4a + - 4<,aSiaS°z*) " K, " «)(™ 



Zq^qlj^qli 



Xi 



u 



(i-m|)(u-m|) 

v a qli a qkj a qki a qlj + a qli a qkj a qki a qlj) ' m Xi m Xj^ 



2( a qU a qkj a qki a qlj + a qh a qkj a qki a qlj) 



(( m l ~ *)( m |* _ *) + ( m li _ ")( m | - «) 



(3.46) 



where 



l,) 2 -4«/2 



4. Numerical results and discussion 



As we assume Xi 2 are lighter than x® 4 and x? is likely to be the LSP, the three types of 
channels: qq — > X1X11 QQ ~~ > X2X2 anc l QQ ~^ X1X2 (associated production) would be the most 
dominant neutralino pair production processes which may lead to the first detection of SUSY 
particles at the LHC||. Here we present the numerical results of the cross sections of these 
three processes and discuss their dependences on the basic SUSY parameters. We divide the 
input MSSM parameters into two parts. One is for for the general parameters included also 
in the SM, and the other is the ino and squark sectors of the MSSM. For the first parameter 



part, we take mz = 91.1887 GeV, sin 2 ^ = 0.2315, a = 1/137.03598. For the second part, 
we just limit the values of Mi, M 2 and /i to be real, positive and below 1 TeV, and take 
tan/5 = 1.5, m fil = mj = 350 GeV, rriu 2 = = 550 GeV, and 6,~ a = = n/A. The three 
ino physical masses are taken as 

m.+ = 150 GeV, m-+ = 550 GeV, m^o = 100 GeV, (4.1a) 

By using Eq.(2.3) with above ino mass values, one may have two choices of parameter 
sets for /i and M 2 , which are in two extreme cases respectively, namely Higgsino-like and 
gaugino-like ino states. For the Higgsino-like case, we get 

M 2 = 534 GeV, \i = 166 GeV, M 1 = 135 GeV, 

and other neutralino masses as 

m^o = 166 GeV, m^o = 184 GeV, m^o = 550 GeV (4.16) 

For the gaugino-like case, we have 

M 2 = 166 GeV, n = 534 GeV, M 1 = 105 GeV, 

and other three neutralino masses as 

m^o = 151 GeV, m^o = 534 GeV, m^o = 554 GeV. (4.1c) 

For complete discussion, we present also the results for the mixture case. We take the 
physical mass inputs as follows for the mixture ino states 

m-+ = 122 GeV, m-+ = 280 GeV, = 101 GeV. 

ii Xi 



the corresponding outputs obtained as 

M 2 = /i = 200 GeV, M x = 150 GeV, 

m^o = 164 GeV, m^o = 201 GeV, m^o = 286 GeV, (4. Id) 

Then the cross sections of the subprocesses qq — ► XiXii XiXlL X2X2 can be numerically eval- 
uated. The cross sections for the three subprocesses are plotted in Fig. 2 (a), (b) and (c) 
respectively, as the functions of the parton-parton collision c.m.s. energy \/§. In Fig. 2(a) 
and (c), all the curves for gaugino-like case show obviously the threshold effects when \/§ 
is just above the threshold energies of the X1X1 an d X2X2 P a i r production, respectively. 
The cross sections for Higgsino-like case are fairly smaller than the corresponding ones for 
gaugino-like case in the process of xlxl P arr production, and the curves for the Higgsino-like 
X2X2 P a i r production are even too small to be plotted in Fig.2(c). On the contrary, the 
Higgsino couplings will enhance abruptly the cross sections for the processes of X1X2 P & i r 
production as shown in Fig.2(b), but the cross sections for gaugino-like and mixture cases 
are negligibly small compared with those for Higgsino-like case. 

The reactions of qq — ► x^Xj are only subprocesses of the parent pp hadron collider. The 
total cross sections of x^Xj P a h productions via qq annihilation in pp collider can be simply 
obtained by folding the cross sections of the subprocesses a[qq — ► XiXj] (l = u,d) with 
the quark and anti-quark luminosity. Adopting quark and anti-quark structure functions 
of the MRS set G given in Ref. || and Q = Vs, we calculate the cross sections a\pp — > 
(R ~ XiXj] (hj — 1) 2) as the functions of the pp collider c.m.s energy yfs, under the same 



10 



conditions of Eq.(4.1). The results are depicted in Fig. 2(d). It is impressive that the X1X2 
pair production rate at the LHC, can even reach as large as 1.5 x 10 2 fb, when the couplings 
are Higgsino-like. We also calculate the X1X2 P a ^ T production at the Tevatron. The cross 
section a[pp — > qq — > X1X2] can De 0-02 pb in Higgsino-like case, when y/s = 2 TeV. We see 
that the light neutralino pair production rates at the Tevatron are far smaller than those at 
the LHC, due to the lower c.m.s energy at Tevatron. 

In Fig. 3(a), (b) and (c), the cross sections a[uu — > XiXi, X1X21 X2X2] versus the parameter 
M 2 , with V§ = 1 TeV and fj, = 400 GeV, are plotted by taking M 1 = 100, 400 GeV and 
1 TeV, respectively. As discussed above with the physical masses m ^+ 2 an d one of m^o, one 
might get three model parameters Mi and M± where M± denote M2 and /i alternatively. 
The properties of neutralinos depend not only on M 2 and ji as charginos do, but also on 
the mass parameter Mi. From Fig. 3, we can say in the case that when there is a large split 
between M± (M + ^> M_), the xlxl P a i r productions via qq collisions have some interesting 
features as follows: 

• When Mi > M±, X\,2 are mainly decide by M 2 and /i: if // M <C Mi M .. Xi,2 
are dominantly composed of Higgsino, and consequently X1X2 P arr production rate is 
significant while XiXi (* — 1)2) pair productions are suppressed; if \x = M_ 3> M 2 = 
M + , Xi is approximately W3 and x 2 is Higgsino, and accordingly X1X1 P & i r production 
is enhanced while X1X2 an< i X2X2 P & i r productions are suppressed. Shown as in Fig. 3 
with the cases of Mi = 1 TeV, and the curve of X2X2 P arr production for the case 
Mi = 1 TeV is too small to be plotted in Fig.3(c). 

11 



• When Mi <C M±, Xi is B. Accordingly, XiXi, X1X2 P & ir productions are dominated 
by gaugino-couplings, and the former is enhance and the latter is suppressed (the 
exceptation for the curve of M\ = 100 GeV in Fig.3(b) are merely because of the mass 
degeneration between m^o^ when M 1 = M 2 = 100 GeV). X2 are decide by M 2 and /i: 
if M 2 = M_ <C n = M + , is mainly W3, and X2X2 P & ir production is enhanced; and 
vice versa. Shown as in Fig. 3 with the cases of Ml = 100 GeV. 

• When M_<Mi<M + , the effects of Mi on Xi 2 are rather weak, and only lead to some 
mixture. Then if jj, = M_ M 2 = M + , Xi,2 are i n some Higgsino-like states, and 
consequently X1X2 P arr production is significant while xI'Xa (i=l,2) pairs are suppressed; 
if M 2 = M_ <C /i = M + , Xi is in gaugino-like state while x\ may be in some mixture 
state, and accordingly x?X? P a i r production is enhanced while X2X2 pairs are 
suppressed. Shown as in Fig. 3 with the cases of Mi = 400 GeV . 

These features can be concluded as that pure gaugino-couplings dominate the production 
of x°iX°i P a ir (as shown by the curves in the areas M 2 < \x = 400 GeV of Fig.3(a,c)), while 
pure Higgsino couplings enhance the xlx'j i} 7^ j) pair productions (as shown by the curve 
of Mi = 1 TeV in the area \x M 2 of Fig.3(b)), and any gaugino-like xi will spoil the large 
production rate of XixX- 

5. Conclusions 

In this work we studied the productions of the light neutralino pairs x?Xi> XiX:L X2X2 y ia 
Drell-Yan process at hadron colliders, and investigated the correlations between the property 
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of the neutralino pair productions and the basic SUSY Lagrangian parameter Mi, M2 and 

From the numerical results, it can be concluded that Higgsino-couplings will increase the 
X1X2 P arr production significantly, while gaugino-couplings enhance the X1X1 an d X2X2 P & i r 
productions. If /i << M 2 and /i^Mi, Higgsino components will be dominated in Xi and 
and the cross section of Higgsino-like X1X2 P & i r production at the LHC, can reach 100 fb. 
Thus, by taking a annual luminosity of pp collision at the LHC being 100 fb' 1 , one can 
accumulate 1 x 10 4 events per year. Therefore the precise measurement of this process is 
suitable in detecting SUSY signals and helpful in determining the basic SUSY parameters. 
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Figure Captions 

Fig.l The Feynman diagrams of the subprocess qq — > XiX^- 

Fig. 2(a) The cross sections of the subprocess qq — > XiXi as functions of v§- The letter 
'u' and 'd' denote the different processes with initial uu and dd collisions respectively, while 
'H', 'g' and 'm' denote the Higgsino-like, gaugino-like and mixture cases respectively. 

Fig. 2(b) The cross sections of the subprocess qq — > X1X2 as functions of vl. The letter 
'u' and 'd' denote the different processes with initial uu and dd collisions respectively, while 
'H', 'g' and 'm' denote the Higgsino-like, gaugino-like and mixture cases respectively. 

Fig. 2(c) The cross sections of the subprocess qq — > X2X2 as functions of VI. The letter 
'u' and 'd' denote the different processes with initial uu and dd collisions respectively, while 
'H', 'g' and 'm' denote the Higgsino-like, gaugino-like and mixture cases respectively. 

Fig. 2(d) The total cross sections of the process pp —>■ qq — ► XiX^j as functions of y/s, 
where 'H', 'g' and 'm' denote the Higgsino-like, gaugino-like and mixture cases, respectively. 
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Fig. 3(a) The cross sections of the subprocess uu 
400 GeV and v 7 ! = 1 TeV . 

Fig. 3(b) The cross sections of the subprocess uu 
400 GeV and = 1 TeV . 

Fig. 3(c) The cross sections of the subprocess uu 
400 GeV and VI = 1 TeV. 



XiXi as functions of M 2 with \i = 



X1X2 as functions of M 2 with = 



X2X2 as functions of M 2 with /i = 
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Fig.2(b) 
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